Photon emission as a source of coherent behaviour of polaritons 
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We show that the combined effect of photon emission and Coulomb interactions may drive an 
exciton-polariton system towards a dynamical coherent state, even without phonon thermalization 
or any other relaxation mechanism. Exact diagonalization results for a finite system (a multilevel 
quantum dot interacting with the lowest energy photon mode of a microcavity) are presented in 
support to this statement. 
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The effectiveness of phonon relaxation in systems of 
excitonic polaritons has been widely discussed recently 
in connection to the possible Bose-Einstein conden- 
sation of polaritons 0- Due to the very small polari- 
ton lifetime, of the order of picoseconds, and the small 
polariton-phonon scattering cross section there is a 
common belief that phonons alone can not account for 
thermalization of the polariton gas to the lattice tem- 
perature (the phonon bottleneck). However, the (time- 
resolved and angle-resolved) observed emission from de- 
caying polaritons Q suggests that, for strong enough 
pumping pulse and a moderate positive detuning, the 
occupation probabilities of polariton states can be fitted 
to a Bose-Einstein distribution with an effective temper- 
ature even lower than the lattice temperature. The ques- 
tion is, thus, what is the source of coherence in such a 
situation? 

In the present paper, we give an answer to the question 
above showing that phonon relaxation is not a necessary 
condition for the polariton system to reach a coherent 
state. Indeed, in our computations we show, for a finite 
system, that the emission of photons is also a source of 
coherence of polaritons leading, under certain conditions, 
to a high occupation probability for the ground state. 

We first consider the model of paper H, in which a 
multilevel quantum dot interacts with the lowest-energy 
photon mode of a microcavity. The Hamiltonian describ- 
ing the system is the following: 
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We include 12 single-electron and 12 single-hole (two- 
dimensional, harmonic oscillator) levels in Eq. We 
assume that both electrons and holes are confined in a 
small spot of a quantum well, in such a way that the 
single-particle spectrum is almost flat: 

T<> e) =E gap , Tj h) = 0. (2) 

The Coulomb coupling constant, /3, is modified at will 
in order to study its effects on the dynamics. (ij\\kl) 
are matrix elements of the Coulomb interaction between 
harmonic oscillator states. The parameter hw gives 
the detuning of the photon energy with respect to the 
(bare) pair energy, and g is the photon-matter coupling 
strength. We consider only states with total angular mo- 
mentum L = 0, which are maximally coupled to the pho- 
ton mode. 

The Hamiltonian, Eq. preserves the polariton 

number, 

N poi = a) a + Y^H h i + e ! e 0A (3) 

i 

and leads to an excitation gap (the gap from the ground 
state to the first excited state of the system with a fixed 
Np i number) roughly proportional to yjN po i g. 

Due to the small, but finite, transparency of the micro- 
cavity mirrors, photons are continuously emitted. In or- 
der to study the evolution of the system due to sucessive 
photon emission events, we first consider the transition 
probabilities: 

Pfi^\^ f \a\^)\ 2 , (4) 
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FIG. 1: Occupation probabilities of the first 20 states in 
the 10-polariton system (initial distribution, panel (a)) and 
in the 5-polariton systems after the emission of 5 photons: 
(b) (3 = 0, hiv = 5 meV, (c) (3 = 2 meV, hco = -3 meV, (d) 
(3 = 2 meV, huj = 5 meV. 
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FIG. 2: The diagonal elements (occupations) as func- 

tions of i (x axis) and N po i (y axis). The areas of the circles 
are proportional to plf pol ' > ■ The left upper panel corresponds 
to the initial distribution, whereas the other three correspond 
to the same sets of parameters as in Fig. The running time 
is such that the mean polariton number is around 6 in each 
case. 



where ipi is a state with polariton number N po i, and ipt 
a state with N po i — 1. The probabilities should be nor- 
malized: 

E p /* = L ( 5 ) 

/ 

We take an initial distribution of weights, p^ 10 '', for the 
states with N po i = 10. After the emission of one photon, 
the weight of the N po i — 9 state -0/ is: 

pf = E^ 10) - (6) 

i 

(8) (7) 

Proceeding in this way, we compute Pf , pf , etc. This 
"dynamics" captures the basics of the Liouville equa- 
tion for the density matrix and, even more, could be 
appropiate for the description of an experiment in which 
the number of photons leaving the cavity is continuously 
measured. 

We draw in Fig. ^the occupation probabilities in the 
5-polariton system after the emission of 5 photons. Panel 



(a) shows the initial distribution in the 10-polariton sys- 
tem. Note that the x-axis represents the 10-polariton 
states, x — 1 refers to the ground state, x — 2 to the 
first excited state, etc. Panels (b) - (d) show the prob- 
abilities of the 5-polariton states for different choices of 
the parameters (3 and Hu>. The photon-matter coupling 
strength, g, was taken equal to 3 meV, corresponding to 
an excitation gap of around 10 meV in the 10-polariton 
system. 

The figure shows that a positive detuning leads to an 
enhancement of the ground-state occupation probability. 
If Coulomb interactions are turned on, the ground state 
occupation reaches 30 % or higher. However, for negative 
detuning the shape of the initial distribution is preserved. 
This fact, coming only from dynamical considerations, 
could be in the basis of the observed polariton coherence 
in Ref. 0. 

These qualitative findings are confirmed by the numer- 
ical solution of the Liouville equation for the density ma- 
trix 0,0: 
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dt 
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3,k 



2^{i>f\<J\i ) j){$j\ a >\'>l ) k)Pki 

2^Pfi(^A^\'^k)(i>k\a\'il>i), 
j,k 



(7) 



where pji — (ipf\p\ipi)- The energy eigenvalues, Ei, and 
transition amplitudes, (ipf\a\4>i), come from the exact 
diagonalization of the Hamiltonian, Eq. (JIJ, for N po i 
ranging from 1 to 10. The parameter n accounts for 
photon losses through the cavity mirrors (Hk w E gap /Q, 
where Q is the cavity quality factor). In our calculations, 
we take k = 0.1 ps -1 . As in the previous example, the 
initial (mixed) state is a distribution of weights, p^f'\ for 
the states of the 10-polariton system. In the simulations, 
we take 20 states (i = 1, . . . 20) in each sector with fixed 
Npoi- 

A sample of the results is shown in Fig. [3 where the 



weights p 



(Npoi) 



for the state i {x axis) and polariton num- 



ber Npoi (y axis) are represented as circles of areas pro- 



portional to p 



(Npoi) 



The distribution of panels in the 



figure is similar to Fig. ^ That is, in the left upper cor- 
ner the initial [t — 0) distribution is drawn. The mean 
number of polaritons is 10 at t = 0. For the other three 
panels, whose (3 and Hlo parameters are the same as in 
Fig. ^ the runing time has been chosen in such a way 
that the mean number of polaritons is 6. 

Besides the same qualitative feature noticed above, 
i.e. enhancement of ground-state occupations when both 
positive detuning and Coulomb interactions are com- 
bined, Fig. [3 reveals that high ground-state occupations 
(as compared with the rest of the states in the same N po i 
sector) is a common characteristic of states with low po- 
lariton numbers, i.e. those who have radiated the most. 

We may define a total ground-state occupation as: 
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From our numerical solution, it follows that p gs 
0.25, 0.11 and 0.07 for the sets of parameters used in Figs. 
□ and i.e. (3 = 2 meV, hw = 5 meV; (3 = 0, ftu = 5 
meV, and (3 = 2 meV, fru; = — 3 meV, respectively. 

It is interesting to draw the time evolution of pg° ■ 
This result is shown in Fig. [3] along with the time de- 
pendence of the mean number of polaritons. The system 
constants are (3 = 2 meV, ftio = 5 meV. At t = 0, we 
have (Np i) — 10, pg t ° ta ^ = 0.036. In around 20 ps, i.e. 
when (Npoi) ~ 8 or two photons have been emitted, the 
total ground-state occupation rises to nearly 0.2. From 
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FIG. 3: The mean number of polaritons (upper panel), and 
the total ground state occupation, defined in Eq. , as func- 
tions of time. The system constants are hu = 5 meV, (3 = 2 
meV. 



t = 20 to 70 ps ((Npoi) from 8 to 5) the total occupa- 
tion still increases, reaching a maximum value of 0.25. 
Afterwards, the total occupation decreases. This later 
decrease is also a dynamical result, not connected with 
the increasing occupation of the vacuum state. Indeed, 
for t = 150 ps, the weight of the vacuum state is only 
0.16. 

In conclusion, we have extracted energy eigenvalues 
and transition amplitudes from exact diagonalization cal- 
culations in order to solve the Liouville equation for the 
density matrix of a decaying polariton system. The re- 
sults show that the emission of photons is an additional 
source of coherence. Further work is needed in order to 
conceptually clarify the role of detuning and Coulomb 
interactions in the dynamical coherence. 

The results should be extended to larger systems in 
order that their relevance to the experiment reported in 
Ref. be more evident. Perhaps, this could be possible 
within the stochastic dynamics framework An inter- 
esting question is whether the total ground-state occupa- 
tion may reach still higher values with increasing number 
of particles. On the other hand, preliminary calculations 
based on a mean-field (BCS) dynamics 9], similar to the 
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semiclassical spin dynamics of the Dicke model [l(j, but 
including Coulomb interactions and photon losses, show 
that the emission of photons through leaky modes of the 
cavity, in spite of its incoherent nature, helps establish- 
ing a dynamical coherent state. Other mechanisms, like 
phonon relaxation, acting at the "horizontal" level in Fig. 
121 could drive the system toward a quasiequilibrium state. 
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